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Abstract 



\Q . We develop a general framework for the construction of probabilities for the time of arrival in 

quantum systems. The time of arrival is identified with the time instant when a transition in the 
detector's degrees of freedom takes place. Thus, its definition is embedded within the larger issue of 
defining probabilities with respect to time for general quantum transitions. The key point in our analysis 
O |. is that we manage to reduce the problem of defining a quantum time observable to a mathematical 

model where time is associated to a transition from a subspace of the Hilbert space of the total system 
to its complementary subspace. This property makes it possible to derive a general expression for the 
probability for the time of transition, valid for any quantum system, with the only requirement that 
the time of transition is correlated with a definite macroscopic record. 

The framework developed here allows for the consideration of any experimental configuration for the 
measurement of the time of arrival, and it also applies to relativistic systems with interactions described 
by quantum field theory. We use the method in order to describe time-of-arrival measurements in high- 
energy particle reactions and for a rigorous derivation of the time-integrated probabilities in particle 
00 ' oscillations. 
t> 

<N 

iri ' 1 Introduction 

O ' 

The time-of-arrival issue [H[2], in its simplest form, amounts to the following problem. One considers an 
initial wave function l^o) for a particle centered around x = and with a positive mean momentum. The 
question is to find the probability P(t)dt that the particle is detected at distance x = L at some moment 
^ between t and t + St. The issue is important in the foundations of quantum mechanics, in relation to 
the role of time in the theory, but also because of the possibility of measuring time-of-arrival probability 
distributions. 

In this paper, we develop a general framework for the construction of probabilities for the time-of- 
arrival in quantum systems. The key idea is the inclusion of the measuring apparatus in the quantum 
description; hence, the time of arrival is defined as a coarse-grained observable associated to the macro- 
scopic records of the apparatus that correspond to a particle's detection. Our framework allows for the 
consideration of any experimental configuration for the measurement of the time of arrival. In particular, 
our method applies to relativistic systems with interactions described by quantum field theory, therefore 
it is particularly suitable for time-of-arrival measurements in high-energy physics. 

A common conclusion emerging from the different approaches to the topic is that the definition of 
probabilities for the time of arrival strongly depends on the specific experimental set-up through which 
the time of arrival is determined. The latter observation is the starting point of our method. Rather 
than attempting to construct a time-of-arrival probability for the properties of the microscopic quantum 
particle, we consider the larger system including the microscopic particle and a measurement apparatus. 
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In the combined system, the time of arrival is associated to a definite macroscopic record of particle 
detection, defined in terms of the apparatus' degrees of freedom. We can then treat the time of arrival 
as a quasi-classical variable, and construct the relevant probabilities using standard expressions from the 
decoherent histories approach to quantum theory [31 SJ EJ EJ [7] . 

We treat the time of arrival as a special instance of the more general notion of the time associated to 
a transition in a quantum system — here the transition refers to the degrees of freedom of the measuring 
apparatus. Thus, in order to describe the time of arrival, we first derive an operator expression for 
the probability associated to transition times in general quantum systems. This probability distribution 
depends only on the Hamiltonian, the projection operators that define the transition and the initial state 
of the system. The derivation requires no assumptions about specific properties of the physical system; 
only that the time of transition is associated to a macroscopic record of observation. Hence, the most 
important contribution of this paper is that it provides a general methodology for the determination of 
the time-of-arrival probabilities associated to an experiment. 

1.1 The time-of-arrival issue 

In quantum theory, probability distributions for observables are constructed using the Born rule, or its gen- 
eralizations. However, the Born rule does not directly apply to the time-of-arrival problem because there 
is no self-adjoint operator for time in quantum mechanics 0[9]. The time t appearing in Schrodinger's 
equation is an external parameter and not an observable. This implies that the squared modulus of the 
time-evolved wave-function \i[)(x,t)\ is not a density with respect to t, and, hence, it cannot serve as a 
definition for the required probabilities [TO] . 

The time-of-arrival problem is a special case of passing from the classical description of a magnitude to 
its associated quantum one. In this case, the methodologies of the classical theory does not generalize to 
the quantum theory. In classical probability theory, time-of-arrival probability distributions are defined in 
terms of probability currents. For a non-relativistic particle of mass to, the probability current associated 
to Schrodinger's equation is 



where \t^t) = e~ lHt \ipo) and H is the free particle's Hamiltonian. However, the current Eq. (TTJ) gives rise 
to negative probabilities for initial states involving superpositions of different momenta [llj . 
The time-of-arrival issue is of significant interest because 

• it provides a simple set-up for exploring different ideas about one of the major foundational issues 
in quantum mechanics, namely, the role of time in the theory. This issue is particularly pertinent on 
the fields of quantum gravity and quantum cosmology: the necessity of a reconciliation between the 
quantum notion of time as an external parameter to a physical system and the dynamical notion 
of time in general relativity generates the so-called problem of time in quantum gravity [12]. In 
particular, ideas about the role of time can be applied both to the time of arrival problem and in 
quantum cosmology HH E] • 

• it is a prototype for other physical problems in quantum theory that involve the definition of 
probabilities with respect to time. Examples include the study of tunneling time [151 116]. that is, of 
the time that it takes a particle to cross a classically forbidden region, and the construction of well- 
defined probabilities for non-exponential decays [17^ I18j. Moreover, the rapid growth of quantum 
information theory brings into the forefront novel issues, such how entanglement is manifested at 
the level of probabilities for time. Conversely, the analogy with photo-detection theory, where 
temporal correlations provide significant information about the electromagnetic-field state [El [2H 
121] , strongly suggests that probabilities and correlations with respect to time can provide novel 
generalized criteria for entanglement. 
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• time-of-arrival probability distributions can in principle be experimentally measured [22] allowing 
for a comparison between different theoretical predictions. The recent OPERA and ICARUS exper- 
iments have determined time-of-arrival probability distribution in the context of neutrino physics 
[23|, 124] . Furthermore, the theory describing particle oscillations |25[ [26] (neutrinos and neutral 
bosons) relies implicitly on the notion of the time of arrival [27] ; the quantity relevant to the exper- 
iments is the total probability of particle detection integrated over all times of arrival t of a particle 
at the detector. 

Currently, there exist several different approaches to the time-of-arrival problem. 

The main limitation of existing approaches is the lack of generality. There is no precise, algorithmic 
procedure allowing for the derivation of the time-of-arrival probabilities for any specific method of particle 
detection. Moreover, investigations are mainly restricted to non-relativistic quantum mechanics. They 
do not incorporate the quantum-field theoretic description of interactions, which would be necessary for 
the study of time-of-arrival in high-energy physics. 

For initial states sharply concentrated in momentum, all approaches lead to probability distributions 
that are peaked around the classical value of the time of arrival. These probability distributions differ in 
their details, the differences being particularly pronounced for initial states with significant momentum 
spread. In fact, time of arrival probabilities are strongly contextual, in the sense that they strongly depend 
on the experimental procedure through which the time of arrival is determined. 

An axiomatic approach, developed by Kijowski, determines the probability distribution P(L, t) for 
the time of arrival of non-relativistic particles by requiring Galilean-covariance and correspondence with 
the classical theory [28]. The resulting expression is 

2 

(2) 

where tpo is the initial state in the momentum representation. For wave functions with support on positive 
values of the momentum p, the probability density Eq. ([2]) is normalized to unity for t G (— oo, oo). Other 
properties of the probability density Eq. ([2]) are discussed in Ref. [29], and generalizations in Ref. [30] . 
Other approaches to the time of arrival include 

• the use of complex potentials modeling the absorption of the particle by a detector at x = L 

[2H S3 EH S3; 

• the consideration of specific detector models [22 [351 ESI EZ] or idealized clocks [35] ; 

• formulations within the decoherent histories approach to quantum mechanics [Tj 140} [39} I4T1 142| [43] ; 

• the analysis the time of first-crossing of x = L for quantum mechanical paths. Such paths are 
defined either using Feynman's prescription [44], or through Bohmian mechanics [15] or through 
phase space quasi-distributions [4"6] . 

1.2 The proposed time-of-arrival algorithm 

In this paper, we address the issue of constructing time-of-arrival probability distributions associated to 
any method of particle detection. To this end, we develop a novel method that is based on the following 
key ideas: (i) the inclusion of the detector degrees of freedom into the quantum description; (ii) the 
definition of the time-of-arrival as a coarse-grained quasi-classical variables associated to macroscopic 
records of particle detection, and (iii) the understanding of the time of arrival as a special case of the 
more general notion of transition-time, that applies to practically all quantum systems. 

Our method is algorithmic: for each experimental set-up, one identifies the operators corresponding to 
the macroscopic records of particle detection, the total Hamiltonian (which includes the Hamiltonian for 
the microscopic particle, the Hamiltonian for the self-dynamics of the detector and an interaction term), 
and the initial state of the combined system. Once the above variables are determined, a unique expression 



P(L,t) 
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for the time-of-arrival probability distribution follows. Moreover, the method involves no restrictions on 
properties of the operators through which the time-of-arrival probabilities are constructed; it applies to 
any quantum system, including relativistic systems interacting through quantum field theory. Hence, it 
is particularly suitable for the study of the time of arrival problem in high-energy physics. 

We do not define the time of arrival as an intrinsic variable characterizing microscopic particles, but 
as a variable associated to degrees of freedom of macroscopic detectors. In particular, we identify the time 
of arrival with the instant t that the apparatus "clicks", that is, with the reading of an external clock 
simultaneous to the creation of a macroscopic record of particle detection. A particle detection constitutes 
a macroscopic, irreversible amplification of a microscopic event; thus, the time-of-arrival t is defined as a 
coarse-grained, quasi- classical variable associated with internal transitions of the detectofi Probabilities 
for such quasi-classical variables are given by the the decoherent histories approach to quantum mechanics. 
We emphasize that in our approach, a particle 's time of arrival is defined only in presence of a definite 
fact of particle detection. 

The time-of-arrival issue then becomes a special case of the broader issue of defining probabilities 
with respect to time for general quantum transitions; other special cases include the definition of decay 
probabilities in unstable systems [T7j , and the construction of coherence functions of the electromagnetic 
field in photo-detection [19|, 121] . The key point in our analysis is that we manage to reduce the 
problem of defining a quantum time observable to a simple mathematical model where time is associated 
to a transition from a subspace of the Hilbert space of the total system to its complementary subspace. 
This property makes it possible to derive a general expression for the probability for the time of transition, 
valid for any quantum system, with the only requirement that the time of transition is correlated with a 
definite macroscopic record. 

Our method involves three steps: (i) the derivation of an expression for the amplitudes associated to 
definite values for the time t that a specific transition takes place in a quantum system; (ii) the construction 
of the associated probabilities through the requirement that the time of transition t is a coarse-grained 
quasi-classical variable; (iii) the specialization to the time-of-arrival problem, by considering a system 
consisting of a microscopic particle interacting with a macroscopic apparatus at distance L from the 
particle source. The result is a general expression for the time-of-arrival probabilities. This expression 
involves certain operators that describe the properties of the particle detector, its interaction with the 
microscopic particle and the type of the recorded observables. These operators are determined by the 
physics of the time-of-arrival experiment in consideration; once they are specified, a unique expression for 
the time-of-arrival probabilities association to the experiment follows. 

Our main results are the following. 

• We construct the time-of-arrival probabilities for three different models of particle detectors. We 
find that, in general, the time of arrival probabilities are strongly dependent on the physics of the 
detector, but there is an important regime where all information about the detector is encoded in 
a single function of the microscopic particle's momentum, the absorption coefficient a(p). 

• We identify a generalization of Kijowski's probability distribution that is valid for any dispersion 
relation for the microscopic particle, including relativistic ones. This probability distribution corre- 
sponds to the case of constant absorption coefficient. 

• We construct a general expression for the time-of-arrival probability in high-energy processes, in 
which the microscopic particle is detected through a reaction described by relativistic quantum field 
theory. In this context, our method leads to a general quantum- measurement-theoretic description 
of particle detectors in high-energy physics. 

• A non-trivial application of our formalism is the rigorous derivation of the time-integrated probabil- 
ities associated to particle-oscillation experiments. We obtain the standard oscillation formulas in 

1 The time-of-arrival variable is coarse-grained in the sense that its value can only be ascertained with macroscopic accuracy, 
and it is quasi-classical in the sense that it is associated to definite macroscopic records. 
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a regime that corresponds to very short values of the decoherence times associated to the particle- 
detection process. Interestingly, we also find that in the regime of larger decoherence times a novel 
non-standard oscillation formula appears. 

The present work originates from ideas in a broader program about the role of time in quantum 
theory [UJ . It employs many concepts from the description of quantum measurements in the decoherent 
histories approach [31 El [6] and it has some similarities to the Davies-Srinivas photo-detection theory [20J . 
Preliminary versions of the method [10] have been employed for the description of tunneling time [16] . 
non-exponential decays [18], and for the study of temporal correlations in particle detectors in relation to 
the Unruh effect [21]. 

The structure of this paper is the following. In Sec. 2, we derive a general formula for the probability 
distribution associated to the time of transition in any quantum system, and then we employ this formula 
in order to define probabilities associated to general time-of-arrival measurements. In Sec. 3, we study 
the time-of-arrival using idealized models for the particle detection process, we obtain a generalization 
of Kijowski's probability distribution, and we apply our formalism to time-of-arrival experiments in high 
energy physics. Sec. 4 contains a non-trivial application of our formalism, in the rigorous derivation of 
the time-integrated probability in particle-oscillation experiments. 

2 Derivation of a general formula for time-of-arrival probabilities 

In this section, we derive an expression for the time-of-arrival probabilities, applicable to any scheme 
of detection for the microscopic particle. To this end, we first construct probabilities associated to the 
time of transition for a general quantum system; the time-of-arrival probabilities arise as a special case. 
In particular, in Sec. 2.1 we present the physical assumptions and the notations of the formalism; in 
Sec. 2.2 we construct the quantum amplitudes associated to a definite value of the time of transition; 
in Sec. 2.3 we construct the probabilities with respect to the time of transition t by assuming that t is 
a coarse-grained quasi-classical variable; and in Sec. 2.4 we consider the special case of time-of-arrival 
measurements where the transitions under consideration correspond to the recording of a particle event 
by a macroscopic detector. 

2.1 Preliminaries 

The key points of our derivation of the probabilities associated to the time of transition of a generic 
quantum system are the following. 

1. Physical transitions in a quantum system are described as transitions between two complementary 
subspaces in the system's Hilbert space. 

2. The time of a transition is defined in terms of a macroscopic record in an apparatus that is correlated 
to the microscopic transition event. Hence, the time of transition is a coarse-grained, decoherent 
observable, for which probabilities can be meaningfully defined. 

Point 1 is standard in ordinary quantum theory. For example, 

• the emission of a photon from an atom corresponds to a transition from the one-dimensional sub- 
space, defined by the electromagnetic field vacuum, to the subspace of single-photon states [20] ; 

• a von Neumann measurement corresponds to a transition from the subspace in which the pointer 
variable X takes its pre-measurement values, to a subspace corresponding to possible measurement 
outcomes; 

• any particle reaction can be described as a transition from the subspace of states associated to the 
initial particles to the subspace of the product particles. 
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Point 2 restricts the context of this study to observed transitions, i.e., to transitions associated to 
definite macroscopic facts of observation. The existence of a macroscopic record necessitates that we 
include a measurement apparatus in the quantum description of the process. 

We proceed to definitions of the quantities relevant to our construction of time-of-transition probabil- 
ities. We denote by % be the Hilbert space of the combined system, describing the degrees of freedom of 
both the quantum system under consideration and the macroscopic apparatus. 

We assume that % splits into two subspaces: H = H + ®H-. The subspace H + describes the accessible 
states of the system given that a specific event is realized; the subspace %- is the complement of H+. For 
example, if the quantum event under consideration is a detection of a particle by a macroscopic apparatus, 
the subspace H + corresponds to all accessible states of the apparatus after a detection event has occurred. 
We denote the projection operator onto 7i + as P and the projector onto 7i- as Q := 1 — P. 

Once the transition has taken place, it is possible to measure the values of variables for the microscopic 
system through their correlation to a pointer variable of the measurement apparatus. We denote by P\ 
projection operators (or, more generally, positive operators) corresponding to different values A of some 
physical magnitude that can be measured only if the quantum event under consideration has occurred. 
For example, when considering transitions associated to the detection of a particle, the projectors P\ may 
be correlated to properties of the microscopic particle, such as position, momentum and energy. The set 
of projectors P\ is exclusive {P\P\> = 0, if A / A'). It is also exhaustive given that the event under 
consideration has occurred; i.e., ^ A P\ = P. 

We also assume that the system is initially (£ = 0) prepared at a state \ipo) £ H+, and that the time 
evolution is governed by the self-adjoint Hamiltonian operator H. 



2.2 Probability amplitudes with respect to the time of transition 

Quantum mechanical probabilities are defined in terms of squared amplitudes. Hence, in order to define 
probabilities for the time of transition, it is necessary to first construct the relevant amplitudes. In 
particular, we employ the definitions of Sec. 2.1, in order to derive the probability amplitude A, [£i, £2]) 
that, given an initial state \ipo), a transition occurs at some instant in the time interval [ii, £2] and a 
recorded value A is obtained for some observable. 

We first consider the case that the relevant time interval is small, i.e., we set t\ = t and £2 = £ + St, 
and we keep only terms of leading to St. Since the transition takes place within the interval [£, £ + St], at 
times prior to £ the state lay within This is taken into account by evolving the initial state \ipo) with 
the restriction of the propagator into that is, with the operator 

S t = lim (Qe~ iAt/N Q) N . (3) 

iV— >-oo 

By assumption, the transition occurs at an instant within the time interval [£, £ + St], after which 
a value A for a macroscopic observable is recorded. This means that in the time-interval [£, £ + St] the 
amplitude transforms under the full unitary operator for time evolution e~ lHSt ~ 1 — iStH. At time £ + St 
the event corresponding to P\ is recorded, so the amplitude is transformed by the action of P\ (or of 

\J~P~\i if P\ is not a projector). For times greater than £ + St, there is no constraint, so the amplitude 
evolves unitarily until some final moment T. 

At the limit of small St, the successive operations above yield 

|Vo; A, [£,£ + ££]> = -iSt e-^ T -^P x HS t \^ ). (4) 

The amplitude \ipol A, [£,£ + St]) is proportional to St. Therefore, it defines a density with respect to 
time: |^ ;A,£) := lim^o ^IV'O - , A, [£, £ + <5£]). From Eq. g]) 

|V>o;A,£) = -ie-^-^PxHSt^o) = -ie~ iAT C(X, £)|V>o>, (5) 
where the class operator C(X,t) is defined as 

C(X,t)=e iflt P x HS t . (6) 
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Since the amplitude |?/>oj A, t) is a density with respect to the time of transition t, its integration with 
respect to t is well-defined. Hence, the total amplitude that the transition occurred at some time in a 
time interval [ii^] is 



W;\,[t 1 ,t 2 ])=-ie~ i6T [ 2 dtC(\,t)\ipo). (7) 

Jtl 

Eq. ([7]) involves the restricted propagator Eq. ([3]), which may be difficult to compute in practice. 
However, there is an important regime where Eq. ([7]) simplifies significantly. We note that if [P, H] = 0, 
i.e., if the Hamiltonian evolution preserves the subspaces T~L±, then \t/jo;X,t) = 0. For a Hamiltonian of 
the form H = Ho + Hj, where [Hq,P] = 0, and Hj a perturbing interaction, we obtain to leading order 
in the perturbation 

C(X,t) = e^ *P A -ff/e- ido *, (8) 
and the restricted propagator Eq. © does not appear in the amplitude Eq. ([7]). 



2.3 Probabilities with respect to the time of transition 



The amplitude Eq. ([5]) squared defines the probability p(X, [ti,t 2 ]) that at some time in the interval [ti, £2] 
a detection with outcome A occurred 



P(A, [t 1} t 2 ]) := (ij};\[hM\^\[hM) 



t2 



t 2 



dt / dt' Tr{e iH ^P x HS\p S t ,HP x ), 



0) 



where p = \ipo){ipo\- 

However, the expression P(X, [ii , ^2]) does not correspond to a well-defined probability measure, be- 
cause it fails to satisfy the Kolmogorov additivity condition for probability measures. To see this, consider 
the probability corresponding to an interval [f 1 , ^3] = [t 1 , ^2] U [^2,^3]- This equals 



P(A, [ti,t 3 ]) = P ( X > h,t 2 ]) + [*2, t 3 ]) + 2Re 



t2 



dt 



t 3 



dt'Tr (c(X,t)p CHX,t') 



t2 



(10) 



Hence, the Kolmogorov additivity condition P(A, [ti , ^3]) = P(X, [^1,^2]) + -P(A, [t 2 , h]), necessary for a 
consistent definition of a probability measure, fails, unless 



2Re 



dt dt'Tr (c{X,t)p C\X,t') 



(11) 



In the consistent/decoherent histories framework, Eq. (jlip is referred to as the consistency condition 
[31 HJ [5] . It is the minimal condition necessary for the association of a consistent probability measure in 
histories. It appears naturally in the present framework, because we construct probabilities associated to 
properties of the system at different moments of time, that is, probabilities associated to histories. 

Eq. (jlip does not hold for generic choices of t\,t 2 and t 3 . However, in a macroscopic system (or in a 
system with a macroscopic component) one expects that Eq. (jlip holds with a good degree of approxi- 
mation, given a sufficient degree of coarse-graining [61 [7J. Thus, if the time of transition is associated to 
macroscopic records in a measurement apparatus, there exists a coarse-graining time-scale a, such that 
the non-additive terms in Eq. (flOl) are strongly suppressed if \t 2 — t\\ » a and | £3 — t 2 \ >> a. Then, 
Eq. ([9]) does define a probability measure when restricted to intervals of size larger than a. 

Hence, assuming a finite coarse-graining time-scale <r, such that Eq. (Ilip is approximately valid for 
I £2 — ill >> o~ an d \t 3 — 1 2 \ >> a, Eq. ([9]) provides a consistent definition of a probability measure for the 
time of transition. 

It is convenient to define the time-of-arrival probabilities by smearing the amplitudes Eq. ([5]) at a 
time-scale of order a rather than employing probabilities for sharply defined time-intervals, as in Eq. ([9]). 
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Then, the time-of-transition probabilities are expressed in terms of densities with respect to a continuous 
time variable. 

To this end, we introduce a family of functions f a (s), localized around s = with width a, and 
normalized so that lim CT _).o fa(s) = S(s). For example, one may employ the Gaussians 

fAs) = vh e ~ A (12) 



The Gaussians Eq. (fT2|) satisfy the following equality. 



s + s' , 



VUt-s)Ut-s') = U{t - -^-)gAs ~ s'), (13) 
where 

g a (s) =exp[- S 2 /(8cx 2 )]. (14) 

Using the functions f a , we define the smeared amplitude \ipo;X,t),j that is localized around the time 
t with width a, as 

\ifo;\,t) ff := f ds^f a (s-t)\i> ;\,s) = f ds^f„(s-t)C(\,s)W ), (15) 
The square amplitudes 



P a {X,t) = (T (^o;A,# Q ;A,i) (7 = J dsdsWf(s - t)f(s> - t)Tr \C(\, s)p C\\, s')\ (16) 
provide a well-defined probability measure: they are of the form Tr[potl(X, t)], where 

fl(A,t) = J dsds'y/Ms - t)f a (s> - t)C\X, s')C(\, s) (17) 



is a density with respect to both variables A and t. 
The positive operator 



oo 



n T (iv) = i - / dt dxn T {x,t), (is) 



o 

corresponds to the alternative N that no detection took place in the time interval [0, oo). n T (iV) together 
with the positive operators Eq. (|17p define a Positive-Operator- Valued Measure (POVM). The POVM 
Eq. (|17h determines the probability density that a transition took place at time t, and that the outcome 
A for the value of an observable has been recorded. 

Using Eq. (fT3|) . and setting S = (s + s')/2, r = s — s', Eq. ([TBI becomes 

P CT (A, t) = f dSUit - S)P(X, t), (19) 

where 

P(X, t) J drg a (r) [d(X, t + ^)po&(\, t - I)] (20) 

Eq. (fT9l) demonstrates that the probability distribution P a is obtained by coarse-graining the classical 
probability distribution P at a scale of a. For systems monitored at a time-scale much larger than a the 
two probability distributions essentially coincide. In that case, the probability density P may be employed 
instead of P a . Moreover, if the resolution scale a is much larger than any timescale characterizing the 
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microscopic system, we can take the limit a — > oo in Eq. (|20p . i.e., we set g a = 1. The resulting probability 
distribution 



P(X, t) = J drTr [c(X, t + ^)po&(\, t - -)] (21) 



is independent of the coarse-graining scale a. 

Eq. (|2ip is a general expression for the time of transition in a quantum systems, which depends only 
on the initial state and the class operators C(A,i), which are constructed solely from the Hamiltonian 
operator and the positive operators associated to the recorded observables. Thus, we have reduced the 
methodological problem of defining time-of-arrival probabilities to a problem of applying a general formula 
to experiments involving time-of-arrival measurements. 

2.4 Derivation of the time-of-arrival probability distribution 

Next, we employ Eq. (|2ip for constructing probabilities associated to time-of-arrival measurements for 
a single particle. To this end, we select the relevant Hilbert space of the theory and the operators that 
appear in Eq. (|2"Tj) . so that they describe the internal transitions of a macroscopic apparatus that take 
place when a microscopic particle is detected. 

The system under consideration consists of a microscopic particle and a macroscopic apparatus. The 
Hilbert space is the tensor product J r ®'H a , where H a is the Hilbert space describing the apparatus's 
degrees of freedom and T is the Fock space 

T = C®Hi@(Hi®Hi) s ,A®---- (22) 

In Eq. (|22p . Hi stands for the Hilbert space associated to a single particle and S and A denote sym- 
metrization and anti-symmetrization respectively. The reason we employ a Fock space is that in many 
detection processes, the microscopic particle is annihilated by the interactions at the detector. Hence, it 
is necessary to consider interactions where the number of microscopic particles is not conserved. 

The Hamiltonian of the total system is H m (g) 1 + 1 (g) H a + Hi n t, where H m describes the dynamics of 
the microscopic particle, H a describes the dynamics of the apparatus and Hint is an interaction term. 

Finally, we specify the macroscopic variables associated to particle detection. These correspond to 
degrees of freedom of the macroscopic apparatus and they are expressed in terms of the positive operators 
1 <8>IIa on T <S>H a , labeled by the values A of a macroscopic observable. The operators Tl\ are defined on 
V. a and they satisfy the completeness relation j d\Il\ = P, where P is the projector onto the subspace 
H+. 

A model for particle detection is defined through the specification of the Hilbert spaces H\ and H m 
and of the operators H m ,H a , Hint and tlx- For each such specification, the POVM Eq. (f2Tj) is uniquely 
constructed and it defines the time-of-arrival probabilities associated to this model. In this sense, the 
method we develop here is fully algorithmic, and it can be employed to construct different models according 
to the different physics of the detection scheme. 

In particular, we assume that the Hamiltonian H m for the microscopic particle is invariant under 
spatial translations. Hence, it depends only on the particle's momentum operator. We denote the particle's 
energy as a function of momentum (the particle's dispersion relation) as e p . We also consider an interaction 
Hamiltonian 



Hint = I <r ' : '' 



x)JUx)+a|(x)Ji(x) , (23) 



where a,(x),dj(x) are the annihilation and creation operators on J s , i labels non-translational degrees 
of freedom (e.g., spin), and ./^(x) are current operators on H a with support in the region D where 
the detector is located. The interaction Hamiltonian Hi n t corresponds to a detection process where the 
microscopic particle is annihilated at the detector. This includes the cases that the particle is absorbed by 
an atom, or that several product particles are created, or a localized excitation is produced. Detection by 
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absorption is not the only possibility, the incoming particle may be detected from energy or momentum 
transfer associated to a scattering process. In this case, the appropriate interaction Hamiltonian is 

H int = J2 I d 3 xA <i (x,x')ot(x')o^(x) 1 (24) 
is 

expressed in terms of a composite operator A^(x, x') defined on the Hilbert space H a . In this paper, we 
shall restrict our considerations to interaction Hamiltonians of the form (|23p . noting that the alternative 
Hamiltonian (I24p pose no particular problem for our approach. 

We place no restriction on the apparatus' Hamiltonian H a , except for the requirement that it does 
not lead to any spurious detection events, i.e. that detection records appear only when particles interact 
with the apparatus. This condition that [H a , P] = 0. It follows that [1 P, H s 1 + 1 H a ] = 0; hence, 
the class operators C(X,t) are obtained from Eq. (jSJ). 

Finally, we consider a single-particle state for the microscopic system 

l^o) = E / d 3 xip 0i (x)a\(x)\0} (25) 

and an initial state of the apparatus Y&o) stationary with respect to the apparatus Hamiltonian H a . 
We assume that the state |^o) satisfies the condition J~(x)\^q) = 0, which guarantees that the only 
transitions in the detector are caused by the interaction with the microscopic particle. We also set the 
scale of energy so that H a \^o) = 0. 

For the class of models specified above, the operator Eq. ([8]) takes the form 



C(X,t) (hfo®>|*o» = Z)/ d 3 xMx,t)\0) ^*V^ii(x)|*o)) , (26) 

where ip(x,t) is the evolution, under e~ tHmt , of the state \tpo) in 7ii in the position representation of a 
single-particle state. 



Then, 



x (y \jr( x yn x e- lH ^'- s ^n x J+(x)\y ). (27) 

In order to simplify the analysis, we ignore the non-translational degrees of freedom of the particle 
and we reduce the system to one spatial dimension, which corresponds to the axis connecting the source 
and the detection region. Substituting Eq. (|27p in Eq. (|2U|) . the probability density takes the form 

P(\,t) = (i> t \S(\)\ip t ) (28) 

where S(X) is an operator defined in terms of its matrix elements in the momentum basis 

x (* 1 J' {x.')\/fhe i6ar \ffhJ + (x) | o ) • (29) 

In a time-of-arrival measurement, the pointer variable A corresponds to position X along the particle's 
axis of motion. We assume the detector is located at a macroscopic distance L from the source, and that 
the accuracy of the detector's position sampling is of order 5. Then, we consider positive operators \/TIl 
corresponding to an unsharp Gaussian sampling of position at X = L, 
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^2 / dXe~ L£ ^\X,a){X,a\, (30) 



The index a in Eq. (|30p refers to the degrees of freedom of the apparatus other than the pointer vari- 
able. Substituting Eq. (130)) into Eq. (j29H . we obtain the general form for the time-of-arrival distribution 
for the particle. 

We note that the familiar probability current Eq. (pQ) for a non-relativistic particle of mass m is of the 
form Eq. (f28|) . with S(L) = \ [p5(x — L) + 5(x — L)p\. However, the operator S(L) of Eq. ([1]) is non-local 
when acting on configuration-space wave- functions, making the probabilities (j28|) non-local functionals of 
the wave- functions tpt( x )- Thus, an interpretation of S(L) as a probability-current operator is untenable. 

Eqs. (I28j) and (|29p define a probability density for the time of arrival that applies to any detection 
scheme, in which the particle is annihilated. The explicit form of the probability density above depends 
on the physics of particle detection. In the following section, we will consider specific models for the 
detector and its interaction with the microscopic particle. 



3 Time-of-arrival probabilities for different models of particle detec- 
tion 

In this section, we explain how Eq. ()28[) is to be employed for the derivation of explicit time-of-arrival 
probabilities associated to specific experiments. We derive such probabilities for different detector models, 
and we demonstrate that there exists a special regime in which all information about the detector is 
encoded in a single function of momentum, the absorption coefficient a{p). From these models, we obtain 
a generalization of Kijowski's probability distribution, valid for any dispersion relation for the microscopic 
particle. Finally, we adapt our formalism for the description of the time-of-arrival in high-energy physics, 
where the detection process involves the creation of several product particles on the detector. 



3.1 Three particle-detection models 

The time-of-arrival probability distribution Eq. (|28j) does not depend only on characteristics of the 
particle (the initial state and the dispersion relation), but also on characteristics of the detector. The 
latter are incorporated into (i) the Hamiltonian H that describes the detector's self-dynamics, (ii) the 
current operator J+(x) that describes the interaction of the microscopic particle with the detector, (iii) 
the initial quantum state of the detector l^o); an d (iv) the set of positive operators 11^ that correspond to 
the pointer variables of the apparatus. Hence, the description of any particular time-of-arrival experiment 
involves a modeling of the detector through a specification of the aforementioned mathematical objects. 
In what follows, we will consider three different types of detector models, and study the properties of the 
associated time-of-arrival probabilities. 



3.1.1 Detection of a coherent particle excitation 

First, we consider the case that the detection of a microscopic particle is accompanied by the creation 
of a particle-like excitation at the detector. The relevant pointer variable is the excitation's position X 
at the locus of the interaction. For example, the excitation may correspond to an excited nucleus or an 
atom. We assume that the pointer-variable X is approximately an autonomous variable, in the sense the 
interaction with the remaining degrees of freedom at the apparatus is weak. Thus, we can identify the 
subspace 7-L a + of H a of accessible states in the apparatus after a particle detection with the Hilbert space 
L 2 (R,dX) of a single particle. 

The Hamiltonian. Assuming negligible interaction with the remaining degrees of freedom, the pointer 
variable evolves unitarily. We consider a Hamiltonian on T-L a + that describes a non-relativistic particle of 
effective mass 
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H a = E + ^-P 2 , (31) 
2/x* 

where P is the excitation's momentum and Eq is a constant corresponding to the energy gap in the 
detector, that is the energy required for the creation of the excitation. 

Initial state and current operator. Since the pointer variable X is correlated to the point x where the 
microscopic particle has been annihilated, we expect that (X\ J + (x)\*S>q) ~ 0, if \x — X\ is significantly 
larger than 5. We, therefore, write {X\ J + (x)\*S?o) = u{x — X), where u(x) is a function that vanishes for 
|x| >> 5. The function u{x) incorporates all information about the initial state of the apparatus and 
its interaction with the microscopic particle. Later, we shall establish that u(x) is closely related to an 
absorption coefficient associated to the detector. 

Pointer variables. For the Hilbert space L 2 (R,dX), the positive operators (|30l) simplify 



1 



dXe —\X){X\. (32) 



(vr^ 2 ) 1 ^ 

With the assumptions above, Eq. (f29l) becomes 

(p'\S(L)\p) = [ dre-^+vW 2 f dxdx' [ dXdX'e' ^'^ ~ L) 



'txS 2 

xu(x - X)u*{x' - X')e ipx - ip ' x 'G(X - X',t), (33) 

where 



G(X-X',r) := (X\e i(Eo+ ^ ) \X') = , /j^ e -» MX 2 ; X ' )2 +^r | (34) 

V — 27TZT 

Carrying out the integrations over x,x' and (X + X')/2 in Eq. (|33p . we obtain 

(p'\S{L)\p) = e- s2 ( p - p ' )2/4 u{p)u*(p')e i(pL - ip ' L 
x / dre-^+vW 2 f dZe-^e i{p+p ' )z/2 G{Z,T), (35) 



where Z = X — X' and u(p) is the Fourier transform of u{x). 

We consider initial states IV'o) emitted from the source at x = with support on positive values of p. 
For macroscopic values of L, the probability density Eq. ([28]) is strongly suppressed if t < 0. Thus, when 
we consider the time-integrated time-of-arrival probability, we can extend the range of integration from 
[0,oo) to (—00,00). We then obtain 



p(l) ■.= r d tp( t ,L) c r dtp(t, L ) = r ^ s , {l } p) \Up)\\ m 

JO J-00 J ^ \Vp\ 

i.e., the time-integrated probability is L-independent. In Eq. ([36]) . is the particle's initial state in the 
momentum representation, and v p = de p /dp is the particle's velocity. 

The time-integrated probability P(L) is a density with respect to L, hence, for a detector of size d in 
the x-direction, the total fraction of detected particles equals P{L)d. For a monochromatic initial state 
with momentum po, |^o| 2 — 2irS(p — po), the total detection probability equals d(po\S(L)\po) /\v Po \. This 
implies that we can define an absorption coefficient a(p) for the detectoJldefined, standardly, as the 
fraction of incoming particles with momentum p > absorbed per unit length of the absorbing medium 

2 The definition of an absorption coefficient depends crucially on the property, established earlier, that the time-integrated 
probability is i-independent. 
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« W = WWW, (37) 

Next, we evaluate the matrix elements of the operator S(L), Eq. f)35f) at the limit Ep^b 2 « 1, that is, 
under the assumption of a very narrow localization of the particle excitation. In this regime, the integral 
over Z in Eq. (|35l) is approximated by 



J dZe-£e^ z G(Z,r) ~ ^Jexp 



<5 2 

— {p + P ') 2 + iE T 



(38) 



Integrate Eq. (|35j) over r, we obtain 



= - n ^^^=u(p)u*(p')e~^ p2+p ' 2) e ipL - ip ' L . (39) 

Eq. (139p simplifies for initial states ipo(p) with mean momentum po such that Eq « e po and with 
momentum spread such that Ap « pq. In this regime, e p + e p i ~ 2^je p e p >. The probability density 
Eq. (gED becomes 



P(L,t) = K 



^e-* 2p2 ^M P y pL -^ 



(40) 



where K = y / Tr5 2 p*/2. 

Eq. ()40p is straightforwardly generalized to initial mixed states of the form 



dpof(po)\ipp ){ippol ( 41 ) 

where \ip po ) denotes an overcomplete set of pure states (for example, coherent states) with mean mo- 
mentum pq and spread Ap, such that Ap « pq, and f(pf0) is a positive- valued function. To see this, 
let as denote the probability density Eq. (|40l) evaluated for initial state \ippo) as P(L,t;ifi p o). Then the 
probability density P{L, t; p) associated to the initial state p, Eq. (|4ip is 

P(L,t;p ) = J dpQf(po)P(L,t;^ p0 ). (42) 

In particular, there is no requirement that the momentum spread of the initial mixed state p be smaller 
than its mean momentum for Eq. (|46p to hold. 

The presence of the cut-off factor e~ &2p2 in Eq. (|40p implies that the detection of particles with 
momenta p » S^ 1 is strongly suppressed. For values of p « S^ 1 , such that the suppression factor 
g -(5 2 p 2 can kg ignored, the current operator S(L) becomes 

S(L) = Ku*{p)hr x l A 8(x - L)h~ 1/A u(p), (43) 

where h is the particle's Hamiltonian. 

The absorbtion coefficient a(p) corresponding to Eq. (1401 is 

*(P) = K^L, (44) 



where v p = de p /dp is the particle's velocity. 
Eq. flHD implies that 

u(p) 



^(ril^y^M (45) 
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where 9(p) is a momentum-dependent phase. The phase 9(p) contributes to the total phase factor i[9(p) + 
pL — e p t] of the integral in Eq. (|40p . The variation of 9{p) is expected to define a microscopic time-scale 
much smaller than L (\9'{p)\ « L), hence, its contribution to the phase of Eq. (|40p is negligible. Then, 
Eq. (pIOl) becomes 

2 

(46) 

Thus, in this regime all information about the measurement apparatus is contained in the absorption 
coefficient a(p). 



P(L,t) 



dp 
2n 



a(p)\v p \4>o(p)e 



ipL- 



3.1.2 Detection of a decoherent excitation 



A more realistic description of detection through the creation of a particle-like excitation involves taking 
into account the influence of the remaining degrees of freedom of the apparatus in the evolution of the 
pointer variable. To this end, we treat the other degrees of freedom as an environment which results to 
non-unitary dynamics on the excitation. The evolution of the pointer variable X is then treated using 
the theory of quantum open systems. 

The model we consider here involves the same Hilbert space, and the same expressions for the initial 
state, the interaction current and the positive operators as the model of Sec. 3.1.1. The difference 
is in the choice of the Hamiltonian. We model the effect of the environment by introducing a stochastic 
term in the Hamiltonian H a of Eq. (|3ip for the pointer variable X, 



where is a Markovian process satisfying 



M[£(i)]=0, 



M[am(t')\ = mt-t'). 



(47) 



(48) 



In Eq. (USD, M. denotes ensemble average and D is a phenomenological constant. 

In this model, the matrix elements of the operator S(L) are of the form Eq. (|33p . where 

G(X - X', t) : = (X'\M[U{s')U^{s)]\X), 



(49) 



where r = s' — s. 

We evaluate the kernel G to second-order in perturbation theory for the noise: 



G(X-X',r) = (X'\e lAT+t ^ T -^ P |r| |X) 



-2iri(r + iD\r\/^ 



_g -2(-r+iD|T|/M«) 



(50) 



The regimes D/n* « 1 of weak coupling to the environment corresponds to leading order in D/fj,* 
to the model studied in Sec. 3.1.1. Here, we consider the opposite limit _D//z* >> 1 of strong coupling to 
the environment. In this regime, the exponential damping behavior dominates in G{X — X',t), and the 
excitation loses rapidly all quantum coherence. In this regime, 



G(X-X',t) ~ 



2D|r| 



-iEqt 



(51) 



Thus, 



dZe 



-JFI+i(P+P')Z/2 



G(Z,t) = (1 + 



D\t\ , 
25^1 



-1/2 



cxp 



5 2 (p + p') 2 



(52) 



In the integration of the right-hand side of Eq. (|52h over r, the dominant contribution comes from values 
of \t\ « Tdec '■= nl$ 2 /D. The time scale Td ec is a decoherence time scale: virtual processes with a 
difference At in the time of arrival, do not contribute to the total probability if At >> Tdec- 
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For sufficiently strong coupling D to the environment, the damping term in Eq. (|52p dominates and 



the right-hand-side term of Eq. ([52]) equals exp 
multiplicative term D ^* p ,y , and Eq. ([33]) becomes 



D(p+p' ) 



r\2 



Integration over r in Eq. ([33]) leads to a 



(P'\S(L)\ P ) = D{ ^ pl? mu%p')e- S2 ^e^ ,L . (53) 

In Eq. (|53|) . contributions from different momenta are suppressed by a factor e~ s2 ( p ~ p '^ . This is the 
most significant qualitative difference from the model of a coherent excitation of Sec. 3.1.1, in which the 
coherence of momentum superpositions is preserved. 

Since the diagonal matrix elements of S{L) are independent of L, the consideration of initial states 
with support on positive values of p, leads to a definition of the absorption coefficient 

«w = i^i^i 2 - < 54 > 

For pure states with mean momentum po and momentum spread Ap, such that Ap << <5 _1 and 
Ap « po, Eq. (|53l) leads again to Eq. (|46l) . with a(p) given by Eq. (i54l) . which is valid for all mixed 
states of the form Eq. (|4ip . 



3.1.3 Particle detection from energy absorption 

The third model to be considered here describes an apparatus where the microscopic particle is detected 
through the excitation of the detector's energy levels. An variation of this model has been studied in detail 
in Ref. |21| . where it was shown that it corresponds to a Glauber-type photo-detector in one regime and 
to a macroscopic Unruh-Dewitt detector in another. The main features of this model are the following. 

Hamiltonian. This model requires no specific form of the Hilbert space % a of the apparatus degrees of 
freedom and of the corresponding Hamiltonian H a . We only assume that the energy eigenstates spanning 
1-L a correspond to a density of states w(E). 

Pointer variable. We consider a single pointer variable X corresponding to the location of the detector, 
and we assume that X is invariant under time translations generated by the Hamiltonian X: [X , H] = 0. 
For example. X may describe different independent sub-detectors, whose location is fixed in space. We 
employ Eq. (|30p for the positive operators corresponding to position sampling. 

Initial state and current operator. Denoting the eigenstates of the Hamiltonian H a by \a), Eq. (|29p for 
the time-of-arrival probabilities involves matrix elements of the form {X, a\ J+(x)\*$>o). We assume that 
these matrix elements depend on the label a only through the energy eigenvalue E(a) associated to a, 
that is we express 

(X,a\J+(x)\V ) =u(x-X,E(a)), (55) 

for some function u(x, E) of position x and energy E. 
With the above assumptions, Eq. (j29l becomes 



(p'\S(L)\p) = V2ire- S2{P i P ' )2 J dEw(E)u 5 (p, E)u* {p 1 , E)e ipL ~ ip ' L 5{E - £p ^ V ). (56) 

The diagonal matrix elements of S(L) are independent of L, leading to the following expression for 
the absorption coefficient 

a(p) = (57) 
\ V P\ 

Again, for pure states with mean momentum p and momentum spread Ap, such that Ap << <5 -1 and 
Ap « p, Eq. ([56]) leads to Eq. (|4"6"]), with a(p) given by Eq. ([57]) . 
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3.2 General expressions for the time-of-arrival probability 

As demonstrated in the three models of Sec. 3.1, the probability densities for the time of arrival depend 
sensitively on the detailed physics of the particle detection scheme. In what follows, we shall show that, 
nonetheless, there exist particular regimes where they take similar form, and that one particular regime 
corresponds to a generalization of Kijowski's probability distribution |28j. 



3.2.1 Initial states with sharp momentum. 



For initial states ipo(p) with mean momentum po and spread Ap and Ap « po, Eqs. (|39j) . (|56p and ([53 
define a probability density of the form 



P(L,t)~K(p c 



dpif>o(p)e 



ipL—ie p t 



+ 0(Ap/p), 



(58) 



where K (po) depends on the special characteristics of each measurement scheme, and it is proportional 
to the absorption coefficient a(po). The total detection probability depends on the measurement scheme, 
however, the time-of-arrival probabilities conditioned upon detection depend only on the initial state. For 
initial states ipo centered around x = 0, the probability density Eq. (f58j) is peaked around the stationary 
phase point L — (de p /dp)t, and corresponds to a mean time of arrival t c i = L/v p . 



3.2.2 The classical regime. 

Before studying the classical limit of the time-of-arrival probabilities derived in Sec. 3.1, we first examine 
a version of the classical time-of-arrival problem. We consider an ensemble of classical free particles, 
described at t = by a phase space probability distribution po(x,p), where x and p are a particle's 
position and momentum, respectively. The time of arrival at x = L is an observable on the state space, 
defined as t = (L — X)/v p . Thus, the probability density for the time of arrival is given by 

P(t) = [ dxdp5(t - -—^) Po (x,p). (59) 



The corresponding equation for quantum mechanical particles is obtained from Eq. (|28j) by expressing 
the initial state ipo m terms of its Wigner function Wo(x,p), 

P(L,t) = J dxdpW s (t,x,p)W (x,p), (60) 

where 

W s (t,x,p) = J |^*(p- ||e^^(L) e -^|p+ |> (61) 

is the Wigner transform of the operator e lHmt S(L)e~ tHmt . 

We compute Eq. (|6ip the operator S(L) of Eq. (|46|) which, as we saw, applies for a large class of initial 
states in different detection schemes. The integral Eq. (f6T|) contains term of the form v p± g , a(p ± |) and 
e , j. Expanding those terms around p and keeping the leading order terms with respect to £, we obtain 

W s (t, X, P) = a(p)S(t - ^^). (62) 

v p 

Hence, the time of arrival distribution becomes 

P(t,L) = [ dxdpa(p)5(t - ^—^)W (x,p). (63) 
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Eq. (|63p corresponds to the classical probability distribution (|59p modified by the absorption coefficient 
a(p) that takes into account the momentum dependence of the particle's detection probability. Eq. ([63]) is 
not normalized to unity; to normalize, one has to divide by the total probability of detection J dtP(t,L). 

We also calculate the first quantum correction by keeping terms proportional to £ 2 in the integral Eq. 
dSU). We obtain 



L — x olq" — (a') 2 + v v v" — (v') 2 „ L — x 
W s (t, X, P) = a(p)S(t -) a 3 PP Kp} d 2 S(t ) + . . . , (64) 

Up oOtUp Dp 



where the primes denote differentiation with respect to the momentum p. 



3.2.3 Generalization of Kijowski's POVM 

The total integral over time of the probability densities for the time of arrival corresponds to the total 
number of particles detected from a detector at x = L. In general, the total probability depends on the 
properties of the initial state. This implies that the normalized probability density P(L,t)/ J dtP(L,t) is 
not a linear functional of the initial state. This is to be expected, since P(L,t)/ J dtP(L,t) corresponds 
to the probability distribution for the time-of-arrival conditioned upon the particle having been detected. 

There is, however, a special regime where the probability density for the time of arrival can be 
normalized to unity by dividing with a constant. This regime corresponds to the domain of validity of 
Eq. (136]) (for all three models of Sec. 3.1), in the special case that the absorption coefficient is a constant. 
In this case, we obtain the probability density 



P(L,t) 



(65) 



which is normalized to unity for all initial states with support on positive values of momentum. 

2 

For a non-relativistic particle (e p = ^— ), Eq. (|65p coincides with Kijowski's POVM for the time-of- 
arrival. Thus, Kijowski's POVM is identified as an ideal time-of-arrival POVM, in which the probability 
of detection does not depend on the incoming particles' momentum. This interpretation is consistent 
with the results of Ref. [10J, where Kijowski's POVM was obtained by modeling the detector as a totally 
absorbing surface for particles of all momenta, defined by Dirichlet boundary conditions at x = L. 

Eq. (|65p provides the generalization of Kijowski's POVM, valid for a general dispersion relation. The 
corresponding current operator is 



S(L) = y/\v\6(x - L)y/\v\, (66) 

where v = (de/dp)(p) is the velocity operator. 

For a relativistic particle of mass m, v = p(p 2 + to 2 ) -1 / 2 ; the probability density Eq. (|65p is not 
invariant under Lorentz transformations. This is to be expected, since the calculation explicitly involves 
the rest-frame of the detector, and the initial state of the detector's degrees of freedom is not Lorentz 
covariant. 



3.3 Time-of-arrival probabilities in high-energy particle reactions 

The simple models studied in Sec. 3.1 allowed us to derive explicit expressions for the time-of-arrival 
probability and led to the generalization of Kijowski's probability distribution. Here, we expand the scope 
of our approach, in order to define time-of-arrival probabilities for experiments in which the particle under 
consideration is detected through high-energy particle reactions. Our aim is to 

1 . Demonstrate that our method applies to high-energy processes and that it can incorporate quantum- 
field theoretic interactions. 

2. Provide a measurement-theoretic characterization of the particle detection process, which, in prin- 
ciple, can lead to explicit mathematical modeling of realistic particle detectors. 
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3. Demonstrate that Eq. ()46p for the time-of-arrival probabilities is of general validity, and not re- 
stricted to the models of Sec. 3.1. 

To this, we consider a measurement scheme where particles, denoted as A are produced from a source 
around x = and propagate towards a detector at distance L from the source, where they are detected 
by means of the process 

A + B 1 + ...B M ^D 1 ... + D N . (67) 

where B m ,D n are particles (different from the A particles), labeled by the indices m = 1,...,M and 
n = 1, . . . , N. The interaction of the A particles with the B m particles produces the particles D n , which 
are the ones that are being detected. Relevant observables are the D n -particles' time of detection, position, 
momentum and so on. These observables are determined through macroscopic pointer variables in the 
detector. 

The Hilbert space Htot associated to the process Eq. (|67|) are described is as a tensor product 
Utot = Ha &> T~L r - Ha is the Fock space J 7 (Hi) corresponding to the ^4-particles. Explicitly, we write 

T(H X ) = c e tfi e (#i ® Hi) s ,a © (Hi ®h x ® Hi) s , a e . . . , (68) 

where S refers to symmetrization (bosons) and A to antisymmetrization (fermions) and Hia is the Hilbert 
space describing a single A particle. The degrees of freedom of the B m and D n particles are incorporated 
into the Hilbert space H r . In general, H r is a tensor product of Fock spaces, one for each field other than 
A, participating in the process Eq. (167p and of a Hilbert space describing other degrees of freedom in the 
detector. 

We first identify the subspaces H± that define the transition under consideration. Since the detection 
proceeds through the process Eq. (|67[) and the A particles are not directly observable, the transition 
corresponds to subspaces of 7i r . The Hilbert space T~L r is decomposed as Hq © Hprod, where Hq is the 
subspace of states prior to the decay A + B m — > D n , and H pro <i is the subspace corresponding to states of 
the decay products. This means that we identify with and H + with H pro d- 

We denote as P the projector P : H. r — > H pro d and Q = 1 — P its complement. The corresponding 
projectors on Utot are 1 ® P and 1 (g) Q. 

Measurements carried out on the product particles correspond to a family of positive operators 1 (g) 
onto Hprod-, where tl\ are positive operators on % r corresponding to different values A of the measured 
observables. They satisfy the completeness relation J d\Tl\ = P. 

The Hamiltonian is of the form 

H = H A ® i + 1 ® H T + H u (69) 

where Ha is the Hamiltonian for free A particles, H r is the Hamiltonian for the B m and D n particles, 
and Hi is the interaction term. 

We assume that any particles B m that are present prior to detection are almost stationary. This 
condition defines our reference frame, that in most cases coincides with the laboratory frame. Assuming 
that the initial state for the B m particles has support to values of momentum much smaller than their 
rest mass, the restriction of the operator H r in the subspace T~Lq is a constant. We choose this constant 
so that the Hamiltonian Ha on single-particle states for the A particles is 

H A = Vm 2 + p 2 -E , (70) 

where p is the A-particle momentum. For simplicity, we ignore the A-particles' spin degrees of freedom. 
The constant Eq is the threshold of the process (|67p 

n m 
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where M^ m and M£> n are the masses of the particles B m and D n respectively 

In the subspace corresponding to the states after the detection events, the Hamiltonian is 



E 



+ f%-M Dn + V[£(t)]), 



(72) 



where p n are the momentum operators for the D n particles. As in the model of Sec. 3.1.2 we treat 
any degrees of freedom of the detector other than the particles involved in the reaction Eq. (167h as an 
environment, and model their action on the product particles by a stochastic term V"[£(i)]. 
The interaction Hamiltonian is 



Hi = Y,J S(x)J+(x) + S+(x)J-(x) 



(73) 



where b,b\ are annihilation and creator operators on Ha, and ./^(x) are current operators, involving 
products of annihilation operators for the B particles and creation operators for the D particles. Since 
no ^4-particles are created during the detection process, the initial state |0o) in H r must satisfy 



0. 



(74) 



With the definitions above, the time-of-arrival probability, restricted to one dimension, is given by 
Eqs. ([28]) and ([29]) . We consider a pointer variable corresponds to the position X\ of the D\ particle. 
Thus, the positive operators IIl of Eq. ([30]) become 



dXt . . .dX n e *^\Xi, . . .,X n )(X!, ...,X n \ 



(75) 



Substituting Eq. ([75]) into Eq. ([29]) . we obtain terms of the form (Xi, . . . X n \ J + (x)\^q) . The positions 
of all particles must be close to the value x of the interaction locus, within the accuracy 5 of Eq. (j32p 
corresponding to the localization of the D\ particle. Hence, we set X\ = X2 = . . . X n and we express 



{X u ... X n \ J+(x)|*o> = S(Xi -X 2 )... 6{Xi - X n )u(X 1 - x), 
in terms of a function u(x). Then, Eq. (I29D becomes 



(76) 



(p'\S(L)\p) 



-<5 2 (p-p') 2 /4 



u(p)u*(p)e 



i(pL—ip'L 



x J dre'^+^f 2 I dZe~^e i{p+p ' )z/2 G{Z,T) 



where 



G(X - X',s' - s) =M 



N 



H(X'\U n (s')Ul(s)\X) 



.n=l 



(77) 



(78) 



In Eq. dZHD, U n (s) = Te (v^+^n+^KWl) stan ds for the evolution operator for the particle D n , and 
M. denotes average over the stochastic process £(•) corresponding to the influence of the environment. 

Eq. (|77p gives the time-of-arrival probabilities associated to a general particle reaction (|67|) . The 
derivation of an explicit formula for a particular detection requires the specification of the field-theoretic 
interaction (as encoded in the function u(p) and a modeling of the term V"[£(i)] that corresponds to the 
stochastic action of the detector's degrees of freedom on the product particles. 

Since the diagonal elements of S(L) in Eq. ()77p are independent of L, we define an absorption 
coefficient of the apparatus as 



a(p) = \u(p)\ 2 / dre 



dZer^e ipZ G(Z,T) 



(79) 
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Following the same arguments and approximations as in Sec. 3.1, we show that for initial states 
narrowly concentrated in momentum, 



P(L,t) 



dp 
2^ 



a(p)\v p \4>o(p)e 



ipL—ie p t 



(80) 



that is, we show that Eq. (|46p for the time-of-arrival probabilities has a broad degree of applicability and 
its derivation requires no special modeling assumptions. 



4 An application: particle oscillations 

We conclude this paper with a non-trivial application of our formalism in the study of particle oscillations. 
We give a rigorous measurement-theoretic derivation of the time-integrated probability associated to 
particle-oscillation experiments. We derive the standard oscillation formula in a regime corresponding to 
very short decoherence time in the detector, and, intriguingly, we identify a different regime (corresponding 
to larger decoherence times) that leads to a novel non-standard oscillation formula. 

Particle oscillations characterize systems in which the single-particle Hilbert space is split into 
subspaces Hi as %\ = ®Hi, such that the dispersion relation e l p is different in each subspace Hi. In 
high-energy physics, particle oscillations appear in neutral bosons and in (massive) neutrinos; the differ- 
ent dispersion relations are due to the different values of the mass mi in the eigenspaces Hi- Particle 
oscillations arise because the creation processes for neutral bosons and neutrinos couple to the particle's 
flavor, and thus produce superpositions of states with different mass. 

To describe particle oscillations, we adapt the formalism of Sec. 3.3. We denote the oscillating particles 
by A, and we consider their detection through the channel f)6T[) that corresponds to the flavor a. The 
Hamiltonian on the single-particle Hilbert space Hi is of the form 



H m = ©; (y p 2 + m 2 - E j , (81) 

where Eq is the threshold of the detection process. 
The interaction Hamiltonian is 

Hi = Y,f d3x pi W U m J+ (x) + SJ (x) U*i J' (x)] , (82) 

i 

where Sj, SJ are annihilation and creator operators on the Fock space Ha, i labels the mass eigenspaces, 
(x) are current operators of flavor a defined on H r , and U a i is the mixing matrhjf). 
The current operator in Eq. (|82p involves products of annihilation operators for the B particles 

and creation operators for the D particles. Since no A-particles are created during the detection process, 

the initial state |<^o) in H r must satisfy 

J-(x)|0 o ) =0. (83) 
We consider a general single-particle state for the A particles 

l^o) =Y, U h I d 3 xbl(^j l0 (x)\0) (84) 
i J 

where |0)a is the vacuum of the Fock space Ha- The dependence the state (f84l on the mixing matrix XJ% 
is a consequence of the creation of the A-particles through an /3-flavor current. 

In particle oscillation experiments, the time of arrival of individual neutrinos is not determined. Thus, 
the relevant quantity is the time-integrated probability density Pp a {L) for the detection along the a-flavor 



3 A more general treatment (relevant to bosons) would involve a kernel instead of a constant for the mixing matrix U a i, 
reflecting the fact that the mixing coefficients may depend on momentum. Here, however, we shall consider initial states 
sharply concentrated in momentum, for which a constant value of U a i provides a good leading-order approximation. 
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channel of ^4-particles created through the /3-flavor channel. Integrating Eq. (|2ip over time t 6 [0, oo), 
we obtain 



/*oo /*oo /» 

p ^ = Ylj Q ds J 

<#o|4(x')V&' f^')^)! V^^WI*o). (85) 



where ^(x, i) is the Schrodinger time evolution of the wave functions ^io( x )) Hz; is given by Eq. ([32]) . 
f7(s) is the evolution operator for the L>-particles of Eq. (|67p including stochastic terms from interaction 
with the environment and Ai denotes stochastic averaging. Reducing the system to one spatial dimension 
and following the procedure that led to the derivation of Eq. ([77]) in Sec. 3.1, we find 

Pp a (L) = UZiUcjUfaUfK J dsj ds' J W^o(p)V5*o(P , )e J(P " P ' )L " i(4s " e; ' S ' ) ^(P,p / , *' - s) (86) 
in terms of the kernel 

/2 
dZe-^e i{p+p ' )z/2 G(Z,T), (87) 

where G(Z, r) is given by Eq. ([78]). In Eq. ([55]). e* = yV + - #o- 
Next, we consider a broad class of initial states 

i/j i0 {x) = foC*)^*, (88) 

where 0o is some real-valued wave function centered around cc = with position spread a x and pi is the 
mean momentum of the state in the i-th mass eigenspace. We assume that \pi — pj\ « \pi\ for all i and 
3- 

In the momentum representation, 

Ao(p) = Mv-Pi), (89) 

where (f>o(j>) is the Fourier transform of 4>o(x). 

In order to perform the integration over p and p' in Eq. (j86H . we expand 

e p = ei+Vi(p- pi), (90) 

where = e^, and V{ = (de i /dp) p= p i . We further assume that the variation of the kernel F(p,p',r) in 
the range of values of momenta where ipm is supported is negligible; hence, F(p,p',r) can be treated as 
a constant in the integration over p and p' in Eq. ()86p . Then, we obtain 

/'OO /*oo 

Pp a {L) = V UyJaU^U^^ / ds / ds'0o(L - Wj -s)^(L - ^ S ')e" i(e " jS " e " lS,) /(s' - a) (91) 
~T Jo Jo 



In Eq. (|9ip . we denoted /(t) = F{p,p,r), where j? is the mean momentum of the state IV'o)- 



The standard oscillation formula. The integrand in Eq. (j9ip is peaked around s = L/v and s' = 

L/vi, namely, around the classical values of the time of arrival corresponding to the dispersion relation in 
the subspace Hi. The probability density Eq. ([9T]) is strongly sensitive on the form of the function f(r), 
which determines whether the amplitudes associated to different times of detection contribute coherently 
in the total probability. In general, the function f(r) depends on the internal dynamics of the detector 
and on the particles' energy scale. It is expected to vanish for sufficiently large values of r. For example, if 
one includes the coarse-graining time-scale a in the derivation, according to Eq (|20p . f(r) would include 
a multiplicative Gaussian term exp(— r 2 /a 2 ) and it would tend to zero for r >> a. In general, the 
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presence of incoherent interactions (as in the model of Sec. 3.1.2) implies that f T is characterized by some 
time-scale Td ec < such that /(r) ~ for r >> t^c- 

The decoherence time scale t^c depends on the physics of the detector and it cannot be specified 
without a precise modeling of the associated interactions, including the effects of the environment.lt plays 
a crucial role to the form of the probability density ([9T]) . If Tdec « \L/vi — L/vj\, amplitudes peaked 
at different values of the time of arrival do not contribute coherently to the total probability. In that 
case, the function f(s' — s) in Eq. (|9ip is effectively proportional to a delta function. It follows that the 
probability density Pp a of Eq. (|9~T1) . is proportional to ds\ip(L, s)\ 2 , as it is assumed in the so-called 
wave-packet description of particle oscillations [HJ 09] . 

In this regime, the evaluation of the probability Eq. (19ip involves the integral 

dsM L ~ Vis)M L ~ v jS )e' i{ ^- ll)s , (92) 

where we extended the limits of integration to (—00,00), since for L >> a x the integrand is strongly 
suppressed for values s < 0. 

The integral Eq. ([92]) is to be estimated subject to the condition \vi — v\ « Vi, for all? and j. To 
this end, we change the integration variable to r = — s + 4, where v is a mean velocity in the initial state 
l^o) 0- The integral Eq. ([9"2]) . then becomes 

e -%-e,)£/« r drM-^L + v ir )M- 6 -^L + v jr )e l ^-^ )r , (93) 
J —00 " " 

where Svi = Vi — v. By Fourier-transforming the functions 4>o in Eq. (|93p . we can estimate the leading- 
order contribution to the integral as <pi - Vj , where 4>i(x) is the inverse Fourier transform of |0(p)| 2 . 
The spread of (j)i(x) is of the same order as a x . Therefore, if 

L « L loc := a x v/ - (94) 

the term <\>\ ^ — - Vj is approximately a constant. The parameter L\ oc is known as the localization length; 

if L > > Lioc the detection probability is strongly suppressed [Ml [50] . 
In the regime where L « Li oc , Eq. (f9Tj) becomes 



Ppa(L) = C^U^U^U^-^^-^, (95) 
u 

where Cfj are positive constants. 

The probability density Eq. ()95[) is a periodic function of the distance L of the detector from the 
source with oscillation wave-numbers 

kji = {pj - pi) - -(ej - e») (96) 

In general, the value of kij in Eq. (196)) depends on the mean values of the momenta pi on the subspaces 
Hi. Since the interactions that produce the oscillations couple to the flavor basis, there is no reason for 
the value of momentum in a given subspace to be consistently larger than the momentum in another 
subspace; hence, we expect that when averaging over the ensemble pi = pj = p. Different assumptions 
about the initial state have been employed in the literature |49} 150] . for example, that the energies are 
equal. In the context of the present formalism, such assumptions do not affect the resulting probability 
distributions. 

Substituting for all pi the mean momentum p of the initial state, we obtain 

2 2 
mf — m~ 

kji = - J ~^i (97) 

and we recover the standard expression for the oscillation wavenumber that applies both to neutral-boson 
and neutrino oscillations. 



4 For example, v may be defined as the arithmetic mean or as the geometric mean of Vi and Vj. As long as \vi — Vj\ << v, 
the results are not affected by the way we choose to define the mean velocity. 
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A non-standard oscillation formula. The standard oscillation formula, Eq. (|97p . was obtained from 
the assumption that « \L/vi — L/vj\. In the opposite regime, where the decoherence time is 
sufficiently large so that T^ ec >> \L/vi — L/vj\, amplitudes peaked at different values of the time of arrival 
contribute coherently in Eq. (I91D . In this regime, the function f(s' — s) is essentially constant and equal 
to /(0). Hence, 

Pp a (L) = £ C^U^U cd U^U^^ < ^ )L , (98) 

ij 

where Cfj are positive constants. In this regime, there is no coherence length L co h- 

The probability density Eq. (198p is a periodic function of the distance L of the detector from the 
source with oscillation wave-numbers 

^fe-ftj-^-^H^-^+^-i-). (s.9) 

Vj Vi Pi Pj \Vi Vj J 

Eq. ([99]) does not depend on the precise choice of the mean momenta pi. Again, setting pi = pj = p, 
Eq. (j99|) becomes 



™ji — _ _ 
p p 



mf + p 2 — J rrij + p- 



(100) 



Eq. (jlOOp is a non-standard oscillation formula. Its dependence on the threshold energy Eq is particularly 
notable, because it implies a different oscillation wavelengths for different channels of detection. 
For Eq = 0, Eq. (fTOOj) becomes 

kij = (m- - rrij)/p, (101) 

i.e., it predicts an oscillation wavenumber twice as large as the standard expression. For neutrinos, this 
non-standard oscillation formula has been derived through other methods — see Ref. [51] and Refs. |491l52j 
for critique. In absence of independent measurements of the mass differences rrii — rrij, the nonstandard 
oscillation formula is indistinguishable experimentally from the standard one. 

For Eq > 0, we consider separately the non-relativistic regime (relevant to neutral bosons) and the 
ultra-relativistic regime relevant to neutrinos. 

In the non-relativistic regime, where \rrii — rrij\ << rrii, we define a "mean mass" m := m '^ mj , to 
obtain 

kji = - - 3 (2m - E ) (102) 

Eq. (I102p has the same dependence on momentum as the standard expression for the oscillation wave- 
length, and thus, it is indistinguishable in absence of an independent measurement of the mass differences 
rrii — rrij . 

In the ultra- relativistic regime, 

('-§)■ ™ 

The momentum dependence of kji in Eq. (I103p differs from that in the standard oscillation formula. The 
difference is more pronounced as p approaches the threshold energy Eq. 

To summarize, the derivation of the standard expression for particle oscillations requires the assump- 
tion that 'interferences' in the time of arrival are decohered at the detector. In absence of a strong 
decoherence effect, the virtual processes, peaked around different values for the time of arrival, contribute 
coherently to the total probability. They result to a different expression of the oscillation wavelength. 
For a given detection process, the standard oscillation formula applies for sufficiently large values of L, 
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and the non-standard oscillation formula applies for sufficiently small values of L. In typical neutrino 
oscillations the baseline L is of the order of 10 2 m and the neutrino energies of the order of hundreds 
of MeVs; then \L/v~i — L/v~j\ is of the order of 10 _22 s. The validity of the standard oscillation formula 
requires that the decoherence time-scale be much smaller than 10~ 22 s. This time-scale is very small, and 
at the moment, there exists no first-principles modeling of actual detectors to establish whether this value 
is physically realistic or not. For this reason, we believe that there is a good prima facie possibility that 
the non-standard oscillation formulas could be physically relevant |53j . 

5 Conclusions 

The main contribution of this paper is the development of an method for determining the time-of-arrival 
probabilities, valid for to any experimental set-up. Our method is algorithmic, in the sense that for 
any modeling of the detector that determines the particle, a unique expression for the time-of-arrival 
probability follows. The method is also general, because it can incorporate any interaction between 
microscopic system and detector, including ones described in terms of relativistic quantum field theory. 
This achieved this result by: (i) reducing the problem of defining quantum temporal observables to a 
mathematical model where time is associated to a transition from a subspace of the Hilbert space of the 
system to its complementary subspace, and (ii) combining a quasi- classical description of the measurement 
records with a fully quantum modeling of the detector's interaction with the microscopic system. 

We constructed time-of-arrival probabilities for three different detector models. We showed that there 
exists a special regime in which all information about the detector is encoded in a single function of 
momentum, the absorption coefficient a{p). From these models, we obtain a generalization of Kijowski's 
probability distribution, valid for any dispersion relation for the microscopic particle. We also adapted 
our formalism for the description of the time-of-arrival in high-energy physics, where the detection pro- 
cess involves the creation of several product particles on the detector. As a non-trivial application of 
the method, we constructed rigorously the time-integrated probability associated to particle-oscillation 
experiments. 

The applicability of the method is not restricted to the time-of-arrival problem. It can be employed 
in order to define probabilities associated to any physical transition (for example, decays of unstable 
systems), provided that such a transition is accompanied by a macroscopic record of observation. Fur- 
thermore, it applies to set-ups involving more than one detectors [21j . Thus, it can be employed for the 
construction of temporal correlation functions between different detectors, or in order to construct tem- 
poral entanglement witnesses associated to measurements in multi-partite systems. Furthermore, since 
the method is compatible with relativistic quantum field theory, it can also be employed towards the 
construction of quantitatively precise models of relativistic quantum measurements [54] . in which the 
spacetime coordinates of an event are random variables rather than externally predetermined parameters. 
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